% begin createhandles.m
% create function handels for OU process
function s = sigmahandle
% Call as:
% s = createhandles
% y1 = s.mycos(1:10);
% y2 = s.mysin(pi);
s.sigma0 = @sigma0;
s.sigma1 = @sigma1;
s.sigma2 = @sigma2;
s.sigma3 = @sigma3;
s.sigma4 = @sigma4;
s.sigma5 = @sigma5;
s.sigma6 = @sigma6;
s.sigma7 = @sigma7;
s.sigma8 = @sigma8;
s.sigma9 = @sigma9;
s.sigma10 = @sigma10;
s.sigma11 = @sigma11;
s.sigma12 = @sigma12;
s.sigma13 = @sigma13;
s.sigma14 = @sigma14;
s.sigma15 = @sigma15;
s.sigma16 = @sigma16;
s.sigma17 = @sigma17;
s.sigma18 = @sigma18;
s.sigma19 = @sigma19;
s.sigma20 = @sigma20;
s.sigma21 = @sigma21;
s.sigma22 = @sigma22;
s.sigma23 = @sigma23;
s.sigma24 = @sigma24;
s.sigma25 = @sigma25;
s.sigma26 = @sigma26;
s.sigma27 = @sigma27;
s.sigma28 = @sigma28;
s.sigma29 = @sigma29;
s.sigma30 = @sigma30;
s.sigma31 = @sigma31;
s.sigma32 = @sigma32;
s.sigma33 = @sigma33;
s.sigma34 = @sigma34;
s.sigma35 = @sigma35;
s.sigma36 = @sigma36;
s.sigma37 = @sigma37;
s.sigma38 = @sigma38;
s.sigma39 = @sigma39;
s.sigma40 = @sigma40;
s.sigma41 = @sigma41;
s.sigma42 = @sigma42;
s.sigma43 = @sigma43;
s.sigma44 = @sigma44;
s.sigma45 = @sigma45;
s.sigma46 = @sigma46;
s.sigma47 = @sigma47;
s.sigma48 = @sigma48;
s.sigma49 = @sigma49;
s.sigma50 = @sigma50;
s.sigma51 = @sigma51;
s.sigma52 = @sigma52;
s.sigma53 = @sigma53;
s.sigma54 = @sigma54;
s.sigma55 = @sigma55;
s.sigma56 = @sigma56;
s.sigma57 = @sigma57;
s.sigma58 = @sigma58;
s.sigma59 = @sigma59;
s.sigma60 = @sigma60;
s.sigma61 = @sigma61;
s.sigma62 = @sigma62;
s.sigma63 = @sigma63;
s.sigma64 = @sigma64;
s.sigma65 = @sigma65;
s.sigma66 = @sigma66;
s.sigma67 = @sigma67;
s.sigma68 = @sigma68;
s.sigma69 = @sigma69;
s.sigma70 = @sigma70;
s.sigma71 = @sigma71;
s.sigma72 = @sigma72;
s.sigma73 = @sigma73;
s.sigma74 = @sigma74;
s.sigma75 = @sigma75;
s.sigma76 = @sigma76;
s.sigma77 = @sigma77;
s.sigma78 = @sigma78;
s.sigma79 = @sigma79;
s.sigma80 = @sigma80;
s.sigma81 = @sigma81;
s.sigma82 = @sigma82;
s.sigma83 = @sigma83;
s.sigma84 = @sigma84;
s.sigma85 = @sigma85;
s.sigma86 = @sigma86;
s.sigma87 = @sigma87;
s.sigma88 = @sigma88;
s.sigma89 = @sigma89;
s.sigma90 = @sigma90;
s.sigma91 = @sigma91;
s.sigma92 = @sigma92;
s.sigma93 = @sigma93;
s.sigma94 = @sigma94;
s.sigma95 = @sigma95;
s.sigma96 = @sigma96;
s.sigma97 = @sigma97;
s.sigma98 = @sigma98;
s.sigma99 = @sigma99;
s.sigma100 = @sigma100;

function y = sigma0(s,T,a,m,sig0)
y = m*(1-exp(-a*s))+exp(-a*s)*sig0;
function y = sigma1(s,T,a,b) 
y = b/a/sqrt(T)*(1-exp(-a*s));
function y = sigma2(s,T,a,b)
y = b*sqrt(2/T)/(a^2+(pi/T)^2)*(a*(cos(pi/T*s)-exp(-a*s))+pi/T*sin(pi/T*s));
function y = sigma3(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(2*pi/T)^2)*(a*(cos(2*pi/T*s)-exp(-a*s))+2*pi/T*sin(2*pi/T*s));
function y = sigma4(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(3*pi/T)^2)*(a*(cos(3*pi/T*s)-exp(-a*s))+3*pi/T*sin(3*pi/T*s));
function y = sigma5(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(4*pi/T)^2)*(a*(cos(4*pi/T*s)-exp(-a*s))+4*pi/T*sin(4*pi/T*s));
function y = sigma6(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(5*pi/T)^2)*(a*(cos(5*pi/T*s)-exp(-a*s))+5*pi/T*sin(5*pi/T*s));
function y = sigma7(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(6*pi/T)^2)*(a*(cos(6*pi/T*s)-exp(-a*s))+6*pi/T*sin(6*pi/T*s));
function y = sigma8(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(7*pi/T)^2)*(a*(cos(7*pi/T*s)-exp(-a*s))+7*pi/T*sin(7*pi/T*s));
function y = sigma9(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(8*pi/T)^2)*(a*(cos(8*pi/T*s)-exp(-a*s))+8*pi/T*sin(8*pi/T*s));
function y = sigma10(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(9*pi/T)^2)*(a*(cos(9*pi/T*s)-exp(-a*s))+9*pi/T*sin(9*pi/T*s));
function y = sigma11(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(10*pi/T)^2)*(a*(cos(10*pi/T*s)-exp(-a*s))+10*pi/T*sin(10*pi/T*s));
function y = sigma12(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(11*pi/T)^2)*(a*(cos(11*pi/T*s)-exp(-a*s))+11*pi/T*sin(11*pi/T*s));
function y = sigma13(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(12*pi/T)^2)*(a*(cos(12*pi/T*s)-exp(-a*s))+12*pi/T*sin(12*pi/T*s));
function y = sigma14(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(13*pi/T)^2)*(a*(cos(13*pi/T*s)-exp(-a*s))+13*pi/T*sin(13*pi/T*s));
function y = sigma15(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(14*pi/T)^2)*(a*(cos(14*pi/T*s)-exp(-a*s))+14*pi/T*sin(14*pi/T*s));
function y = sigma16(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(15*pi/T)^2)*(a*(cos(15*pi/T*s)-exp(-a*s))+15*pi/T*sin(15*pi/T*s));
function y = sigma17(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(16*pi/T)^2)*(a*(cos(16*pi/T*s)-exp(-a*s))+16*pi/T*sin(16*pi/T*s));
function y = sigma18(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(17*pi/T)^2)*(a*(cos(17*pi/T*s)-exp(-a*s))+17*pi/T*sin(17*pi/T*s));
function y = sigma19(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(18*pi/T)^2)*(a*(cos(18*pi/T*s)-exp(-a*s))+18*pi/T*sin(18*pi/T*s));
function y = sigma20(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(19*pi/T)^2)*(a*(cos(19*pi/T*s)-exp(-a*s))+19*pi/T*sin(19*pi/T*s));
function y = sigma21(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(20*pi/T)^2)*(a*(cos(20*pi/T*s)-exp(-a*s))+20*pi/T*sin(20*pi/T*s));
function y = sigma22(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(21*pi/T)^2)*(a*(cos(21*pi/T*s)-exp(-a*s))+21*pi/T*sin(21*pi/T*s));
function y = sigma23(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(22*pi/T)^2)*(a*(cos(22*pi/T*s)-exp(-a*s))+22*pi/T*sin(22*pi/T*s));
function y = sigma24(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(23*pi/T)^2)*(a*(cos(23*pi/T*s)-exp(-a*s))+23*pi/T*sin(23*pi/T*s));
function y = sigma25(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(24*pi/T)^2)*(a*(cos(24*pi/T*s)-exp(-a*s))+24*pi/T*sin(24*pi/T*s));
function y = sigma26(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(25*pi/T)^2)*(a*(cos(25*pi/T*s)-exp(-a*s))+25*pi/T*sin(25*pi/T*s));
function y = sigma27(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(26*pi/T)^2)*(a*(cos(26*pi/T*s)-exp(-a*s))+26*pi/T*sin(26*pi/T*s));
function y = sigma28(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(27*pi/T)^2)*(a*(cos(27*pi/T*s)-exp(-a*s))+27*pi/T*sin(27*pi/T*s));
function y = sigma29(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(28*pi/T)^2)*(a*(cos(28*pi/T*s)-exp(-a*s))+28*pi/T*sin(28*pi/T*s));
function y = sigma30(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(29*pi/T)^2)*(a*(cos(29*pi/T*s)-exp(-a*s))+29*pi/T*sin(29*pi/T*s));
function y = sigma31(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(30*pi/T)^2)*(a*(cos(30*pi/T*s)-exp(-a*s))+30*pi/T*sin(30*pi/T*s));
function y = sigma32(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(31*pi/T)^2)*(a*(cos(31*pi/T*s)-exp(-a*s))+31*pi/T*sin(31*pi/T*s));
function y = sigma33(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(32*pi/T)^2)*(a*(cos(32*pi/T*s)-exp(-a*s))+32*pi/T*sin(32*pi/T*s));
function y = sigma34(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(33*pi/T)^2)*(a*(cos(33*pi/T*s)-exp(-a*s))+33*pi/T*sin(33*pi/T*s));
function y = sigma35(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(34*pi/T)^2)*(a*(cos(34*pi/T*s)-exp(-a*s))+34*pi/T*sin(34*pi/T*s));
function y = sigma36(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(35*pi/T)^2)*(a*(cos(35*pi/T*s)-exp(-a*s))+35*pi/T*sin(35*pi/T*s));
function y = sigma37(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(36*pi/T)^2)*(a*(cos(36*pi/T*s)-exp(-a*s))+36*pi/T*sin(36*pi/T*s));
function y = sigma38(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(37*pi/T)^2)*(a*(cos(37*pi/T*s)-exp(-a*s))+37*pi/T*sin(37*pi/T*s));
function y = sigma39(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(38*pi/T)^2)*(a*(cos(38*pi/T*s)-exp(-a*s))+38*pi/T*sin(38*pi/T*s));
function y = sigma40(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(39*pi/T)^2)*(a*(cos(39*pi/T*s)-exp(-a*s))+39*pi/T*sin(39*pi/T*s));
function y = sigma41(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(40*pi/T)^2)*(a*(cos(40*pi/T*s)-exp(-a*s))+40*pi/T*sin(40*pi/T*s));
function y = sigma42(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(41*pi/T)^2)*(a*(cos(41*pi/T*s)-exp(-a*s))+41*pi/T*sin(41*pi/T*s));
function y = sigma43(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(42*pi/T)^2)*(a*(cos(42*pi/T*s)-exp(-a*s))+42*pi/T*sin(42*pi/T*s));
function y = sigma44(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(43*pi/T)^2)*(a*(cos(43*pi/T*s)-exp(-a*s))+43*pi/T*sin(43*pi/T*s));
function y = sigma45(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(44*pi/T)^2)*(a*(cos(44*pi/T*s)-exp(-a*s))+44*pi/T*sin(44*pi/T*s));
function y = sigma46(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(45*pi/T)^2)*(a*(cos(45*pi/T*s)-exp(-a*s))+45*pi/T*sin(45*pi/T*s));
function y = sigma47(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(46*pi/T)^2)*(a*(cos(46*pi/T*s)-exp(-a*s))+46*pi/T*sin(46*pi/T*s));
function y = sigma48(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(47*pi/T)^2)*(a*(cos(47*pi/T*s)-exp(-a*s))+47*pi/T*sin(47*pi/T*s));
function y = sigma49(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(48*pi/T)^2)*(a*(cos(48*pi/T*s)-exp(-a*s))+48*pi/T*sin(48*pi/T*s));
function y = sigma50(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(49*pi/T)^2)*(a*(cos(49*pi/T*s)-exp(-a*s))+49*pi/T*sin(49*pi/T*s));
function y = sigma51(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(50*pi/T)^2)*(a*(cos(50*pi/T*s)-exp(-a*s))+50*pi/T*sin(50*pi/T*s));
function y = sigma52(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(51*pi/T)^2)*(a*(cos(51*pi/T*s)-exp(-a*s))+51*pi/T*sin(51*pi/T*s));
function y = sigma53(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(52*pi/T)^2)*(a*(cos(52*pi/T*s)-exp(-a*s))+52*pi/T*sin(52*pi/T*s));
function y = sigma54(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(53*pi/T)^2)*(a*(cos(53*pi/T*s)-exp(-a*s))+53*pi/T*sin(53*pi/T*s));
function y = sigma55(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(54*pi/T)^2)*(a*(cos(54*pi/T*s)-exp(-a*s))+54*pi/T*sin(54*pi/T*s));
function y = sigma56(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(55*pi/T)^2)*(a*(cos(55*pi/T*s)-exp(-a*s))+55*pi/T*sin(55*pi/T*s));
function y = sigma57(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(56*pi/T)^2)*(a*(cos(56*pi/T*s)-exp(-a*s))+56*pi/T*sin(56*pi/T*s));
function y = sigma58(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(57*pi/T)^2)*(a*(cos(57*pi/T*s)-exp(-a*s))+57*pi/T*sin(57*pi/T*s));
function y = sigma59(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(58*pi/T)^2)*(a*(cos(58*pi/T*s)-exp(-a*s))+58*pi/T*sin(58*pi/T*s));
function y = sigma60(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(59*pi/T)^2)*(a*(cos(59*pi/T*s)-exp(-a*s))+59*pi/T*sin(59*pi/T*s));
function y = sigma61(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(60*pi/T)^2)*(a*(cos(60*pi/T*s)-exp(-a*s))+60*pi/T*sin(60*pi/T*s));
function y = sigma62(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(61*pi/T)^2)*(a*(cos(61*pi/T*s)-exp(-a*s))+61*pi/T*sin(61*pi/T*s));
function y = sigma63(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(62*pi/T)^2)*(a*(cos(62*pi/T*s)-exp(-a*s))+62*pi/T*sin(62*pi/T*s));
function y = sigma64(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(63*pi/T)^2)*(a*(cos(63*pi/T*s)-exp(-a*s))+63*pi/T*sin(63*pi/T*s));
function y = sigma65(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(64*pi/T)^2)*(a*(cos(64*pi/T*s)-exp(-a*s))+64*pi/T*sin(64*pi/T*s));
function y = sigma66(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(65*pi/T)^2)*(a*(cos(65*pi/T*s)-exp(-a*s))+65*pi/T*sin(65*pi/T*s));
function y = sigma67(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(66*pi/T)^2)*(a*(cos(66*pi/T*s)-exp(-a*s))+66*pi/T*sin(66*pi/T*s));
function y = sigma68(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(67*pi/T)^2)*(a*(cos(67*pi/T*s)-exp(-a*s))+67*pi/T*sin(67*pi/T*s));
function y = sigma69(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(68*pi/T)^2)*(a*(cos(68*pi/T*s)-exp(-a*s))+68*pi/T*sin(68*pi/T*s));
function y = sigma70(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(69*pi/T)^2)*(a*(cos(69*pi/T*s)-exp(-a*s))+69*pi/T*sin(69*pi/T*s));
function y = sigma71(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(70*pi/T)^2)*(a*(cos(70*pi/T*s)-exp(-a*s))+70*pi/T*sin(70*pi/T*s));
function y = sigma72(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(71*pi/T)^2)*(a*(cos(71*pi/T*s)-exp(-a*s))+71*pi/T*sin(71*pi/T*s));
function y = sigma73(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(72*pi/T)^2)*(a*(cos(72*pi/T*s)-exp(-a*s))+72*pi/T*sin(72*pi/T*s));
function y = sigma74(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(73*pi/T)^2)*(a*(cos(73*pi/T*s)-exp(-a*s))+73*pi/T*sin(73*pi/T*s));
function y = sigma75(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(74*pi/T)^2)*(a*(cos(74*pi/T*s)-exp(-a*s))+74*pi/T*sin(74*pi/T*s));
function y = sigma76(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(75*pi/T)^2)*(a*(cos(75*pi/T*s)-exp(-a*s))+75*pi/T*sin(75*pi/T*s));
function y = sigma77(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(76*pi/T)^2)*(a*(cos(76*pi/T*s)-exp(-a*s))+76*pi/T*sin(76*pi/T*s));
function y = sigma78(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(77*pi/T)^2)*(a*(cos(77*pi/T*s)-exp(-a*s))+77*pi/T*sin(77*pi/T*s));
function y = sigma79(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(78*pi/T)^2)*(a*(cos(78*pi/T*s)-exp(-a*s))+78*pi/T*sin(78*pi/T*s));
function y = sigma80(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(79*pi/T)^2)*(a*(cos(79*pi/T*s)-exp(-a*s))+79*pi/T*sin(79*pi/T*s));
function y = sigma81(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(80*pi/T)^2)*(a*(cos(80*pi/T*s)-exp(-a*s))+80*pi/T*sin(80*pi/T*s));
function y = sigma82(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(81*pi/T)^2)*(a*(cos(81*pi/T*s)-exp(-a*s))+81*pi/T*sin(81*pi/T*s));
function y = sigma83(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(82*pi/T)^2)*(a*(cos(82*pi/T*s)-exp(-a*s))+82*pi/T*sin(82*pi/T*s));
function y = sigma84(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(83*pi/T)^2)*(a*(cos(83*pi/T*s)-exp(-a*s))+83*pi/T*sin(83*pi/T*s));
function y = sigma85(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(84*pi/T)^2)*(a*(cos(84*pi/T*s)-exp(-a*s))+84*pi/T*sin(84*pi/T*s));
function y = sigma86(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(85*pi/T)^2)*(a*(cos(85*pi/T*s)-exp(-a*s))+85*pi/T*sin(85*pi/T*s));
function y = sigma87(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(86*pi/T)^2)*(a*(cos(86*pi/T*s)-exp(-a*s))+86*pi/T*sin(86*pi/T*s));
function y = sigma88(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(87*pi/T)^2)*(a*(cos(87*pi/T*s)-exp(-a*s))+87*pi/T*sin(87*pi/T*s));
function y = sigma89(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(88*pi/T)^2)*(a*(cos(88*pi/T*s)-exp(-a*s))+88*pi/T*sin(88*pi/T*s));
function y = sigma90(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(89*pi/T)^2)*(a*(cos(89*pi/T*s)-exp(-a*s))+89*pi/T*sin(89*pi/T*s));
function y = sigma91(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(90*pi/T)^2)*(a*(cos(90*pi/T*s)-exp(-a*s))+90*pi/T*sin(90*pi/T*s));
function y = sigma92(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(91*pi/T)^2)*(a*(cos(91*pi/T*s)-exp(-a*s))+91*pi/T*sin(91*pi/T*s));
function y = sigma93(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(92*pi/T)^2)*(a*(cos(92*pi/T*s)-exp(-a*s))+92*pi/T*sin(92*pi/T*s));
function y = sigma94(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(93*pi/T)^2)*(a*(cos(93*pi/T*s)-exp(-a*s))+93*pi/T*sin(93*pi/T*s));
function y = sigma95(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(94*pi/T)^2)*(a*(cos(94*pi/T*s)-exp(-a*s))+94*pi/T*sin(94*pi/T*s));
function y = sigma96(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(95*pi/T)^2)*(a*(cos(95*pi/T*s)-exp(-a*s))+95*pi/T*sin(95*pi/T*s));
function y = sigma97(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(96*pi/T)^2)*(a*(cos(96*pi/T*s)-exp(-a*s))+96*pi/T*sin(96*pi/T*s));
function y = sigma98(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(97*pi/T)^2)*(a*(cos(97*pi/T*s)-exp(-a*s))+97*pi/T*sin(97*pi/T*s));
function y = sigma99(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(98*pi/T)^2)*(a*(cos(98*pi/T*s)-exp(-a*s))+98*pi/T*sin(98*pi/T*s));
function y = sigma100(s,T,a,b) 
y = b*sqrt(2/T)/(a^2+(99*pi/T)^2)*(a*(cos(99*pi/T*s)-exp(-a*s))+99*pi/T*sin(99*pi/T*s));
% end createhandles.m